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Abstract 

We develop a general approach to calculations of the electron spectrum of metal- 
lic multi-wall carbon nanotubes (MWNT) with arbitrary number of coaxial layers. 
It is based on the model with singular attractive potential of equidistant conductive 
cylinders. The knowledge of one-electron spectrum allows to construct the corre- 
sponding Green function and then to calculate the entropy and density of states 
for MWNT. We analyze the tunnelling between the nanotube and normal metal 
electrode. The possibility of direct determination of one-electron density of states 
by measurements of the tunnelling conductivity at low temperatures is proved and 
the necessary restrictions on temperature are formulated. We discuss briefly the 
conflicting experimental observations of electronic properties of MWNT. 

PACS: 73.22.f; 76.63.Fg. 
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1 Introduction 

Recently unusual properties of carbon nanotubes attract much interest. These mesoscopic 
systems have demonstrated a remarkable interplay of dimensionality, interaction and dis- 
order pp. The well-pronounced one- dimensional character of electron motion in single- wall 
carbon nanotubes (SWNT) with strongly separated levels of angular quantization led to 
attempts to describe their behavior in terms of the Luttinger liquid model [3 EH HI E] • 
As soon as metallic SWNT has two conducting channels, the ID Luttinger liquid model 
may apply to this case when the interchannel coupling is negligible. 

Experimentally, several contradictory observations of the electronic properties of SWNT 
at low temperatures were made. Such an evidence of the Luttinger liquid behavior as a 



power law divergence of resistance when temperature fell down to 10 K^, was recently 
attributed to Coulomb blockade in tunnel junctions [B]. The measured with the use of 
low-resistive metallic contact low-temperature resistance of individual SWNT, [Hj does 
not manifest traces of Luttinger liquid behavior. 

The multiwall nanotubes (MWNT), composed of several concentric graphite shells 
show properties which are clearly consistent with the weak-localization features of diffusive 
transport in magnetoconductivity and zero-bias anomaly in the tunnelling density of 
states [7J |H1 El- The intertube coupling in MWNT was demonstrated to be essential for 
interpretation of experimental data on electric transport properties of MWNT [TT1 ITU] . 
On the other hand, if electrods contact only the outermost tube in MWNT, it behaves 
as an individual cylinder with the diameter an order of magnitude larger than that of 
SWNT and may manifest ballistic properties |12j . 

In the present paper we study electron states in MWNT. In this case, contrary to 
SWNT, the inter-tube electron hopping eliminates the specific one-dimensional proper- 
ties and we may solve the corresponding spectral problem within the standard quantum 
mechanical many-electron approach. In order to describe the electron spectrum of MWNT 
one has to take into consideration that the motion of electron around the circumference of 
individual tube is quantized. Due to the periodic boundary conditions an integer number 
of wavelengths must fit around the tube. Along the tube, however, electronic states are 
not restricted and electrons can move ballistically. Because of the circumferential quan- 
tization, the electron states in the tube do not form a single wide energy band. Instead, 
they split into a number of one-dimensional sub-bands, with band onsets at different ener- 
gies. In MWNT an additional effect of intertube interaction appears. It results in further 
splitting of energy levels. The characteristic distance between the levels is governed by 
the probability of electron tunnelling between adjacent tubes. 

We model the metallic MWNT as a set of equidistant coaxial cylinders with attractive 
Kronig-Penney type potential along the radius and study the electron states of such a 
system. We develop an universal formalism for numerical evaluation of electron energy 
levels. It allows to calculate energy level splitting due to intertube hopping for MWNT 
with any number of coaxial tubes. Then we consider the filling of the obtained band 
structure by electrons and calculate such electronic properties as the one-particle density 
of states, electronic contribution to entropy and, finally, the tunnelling conductivity of 
MWNT. We also prove the possibility of direct determination of the electron spectrum in 
MWNT by means of the tunnelling microscopy and formulate the corresponding criterion 
on experimental temperatures. 

2 Electrons on coaxial nanotubes 

2.1 The model and notation 

We shall consider the system of graphite nanotubes as a set of equidistant coaxial con- 
ducting cylinders. The number of cylinders iV is arbitrary. Electrons are held to the 
cylinders by a singular attractive potential. The single electron Hamiltonian takes the 
form: 

h 2 v 2 N 

H = \ U5(r — ri), with r» = r\ + (i — l)a. (1) 
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Here N is the number of coaxial tubes and r, are their radii. The radius of the inner tube 
is ri and the adjacent tubes are separated by distance a < r±. We shall assume that the 
potential is strong. Quantitatively this means that U ^> -^-y. 
The Schrodinger equation for electron reads: 

N 

y-J2 US ( r - r i)V = E tf . (2) 

8=1 

We shall solve the problem in cylindrical coordinates with z-axis pointing along the tube. 
Let us take the electron wave function in the form 

V(z, r, <j>) = y^ i/j n (r) exp ( -p z z + in§ J , (3) 

n=0 ^ ' 

The electron has momentum p 2 and moves freely along the tube. The angular motion of 
the electron is characterized by magnetic quantum number n. The potential affects only 
the radial part of the wave function: 

dr 2 r dr + r 2^ h2 l^o\r r t )ip n ii<, 2m )Vn. W 

Note that the electrons are localized on the tubes and the right hand side of the equation 
is negative. 

The potential energy may be characterized by the effective width of 5-wells: 

\/mU 

It may be seen from the corresponding one- dimensional problem that A controls the spread 
of wave functions out of the tubes. Large values of U correspond to the case A <C a, r. 

Another important quantity is the spectral parameter. It has the dimension of inverse 
length: 

K = - \Jp 2 z — 2mE. (6) 

I v 

The spectral parameter sets the radial scale for wave function and its values determine 
the discrete component of energy spectrum. 

Our goal is to find k as a function of A depending on the system geometry. After the 
standard rescaling x = nr the equation for ip n (x) takes the form: 

2 N 
x 2 i\)" n + xip' — (x 2 + n 2 )il) n = -y x 2 5(x — xi) ip n , with Xi — KTi. (7) 

i=l 

We shall describe square integrable solutions of this equation in the next section. 

2.2 Matching conditions 

The potential on the right hand side of equation is zero almost everywhere excluding 
the points Xi = KTi. Between the tubes the equation is of the modified Bessel type. Hence 



on the intervals x« < x < Xj+i the wave function ip n must be a superposition of modified 
Bessel functions: 

ip n (x) = Afl n {x) + B?K n {x), for x e [a:*, x m ], x = 0. (8) 

The 2N + 2 coefficients A™, Bf, i = 0, . . . N, are determined by matching conditions at 
singular points of potential (0 ITTjl . boundary conditions (fT2"|) and overall normalization 
of the wave function (|T5|) . 

The first iV relations between A™_ t , B™_ x and A", Bf on the next interval come from 
the continuity of solutions, ij) n {xi — 0) = i^ n (xi + 0). Substituting (jHJ) for ip n we get: 

AtxWi) + BZ_ x K n { Xi ) = A?I n ( Xi ) + B?K n (x t ), i = l,...N. (9) 

The second set of conditions is obtained from matching the derivatives ip' n . They 
change at singular points Xi stepwise and if/^ behaves like (^-function. Integrating the two 
sides of equation (jJJ) over an infinitesimal interval (x» — e, Xi + e) we obtain for e — > 0: 



lim / x 2 < dx = x 2 Vnl*!-o = r^n{xi) and 






Xi+0 r, 



Xi-0 



kA 



(10) 



This leads to iV more equations of the form: 

{A? - A?_M Xi ) + (B? - Bf^K'Jxi) _ 2 
AfUxi) + B?KJxi) k\ 



1,...N, (11) 



where we have used (JHJ) in the denominator. 

The requirement of regularity of the wave function at zero and infinity results in two 
more relations. Remember that the modified Bessel function I„(x) is regular at x = 
and grows at x — > oo whereas K n (x), on the contrary, is singular at x = but goes to 
zero at infinity. Hence one must set 

A% = B- = 0. (12) 

It may be easily seen that the listed 2N + 2 equations are not independent. Namely, 
rescaling all A — > aA, B —>■ aB does not break the equalities. The missing requirement 
that helps to fix the values of A and B is the normalization condition 



\ijj n (x)\ xdx— \ijj n (r)\ r dr = 1. (13) 

Jo 

We shall show how to solve equations (0 El E} i n the next section. 

2.3 The recursions 

The standard way to find the spectrum of the described system of homogenous equations 
would be to study the dependence of its determinant on spectral parameter. We propose 
a recursive procedure that allows to deduce the spectral equation without calculating the 
determinant. 

The wave functions of electrons in the system of iV coaxial nanotubes are given by 
formula (JHJ) where 2N + 2 coefficients A™ and B™ may be found up to a factor from the 



following system of 2N + 2 homogenous linear equations (we summarize once more (0 

mm): 

(X? - AU)In{xi) - (3f-i - B?)K n {zi) = 0, (14a) 

i = 1,...N; 

(A«-^ 1 K(^ + (i?r-^ 1 )K(^ + A (A nj n(Xi)+i?r ^ (Xj)) = 0, (14b) 

i = 1,...N; 

Bl = 0; (14c) 

A n N = (14d) 



This system has nonzero solutions only provided that its rank is less than 2N + 2. 
Thus we arrive at a symbolic equation det||LHS|| = 0. The determinant of the LHS 
is an iV-th order polynomial of (A/t) _1 with complicated (due to the presence of Bessel 
functions) coefficients that depend on k and on geometrical parameters of the system. 
Zeros of this function with respect to k determine the energy spectrum of the system. 

An alternative way to obtain the spectral equation without calculating the determi- 
nant is to eliminate variables. This may be done recursively. Equations ()14al I14b|) link 
coefficients Af_ ± , B™_ x to A™, B™ on the next interval. Applying them successively one 
can express all A™, B™ in terms of Aq, Bq. At the end, after imposing the boundary 
conditions (|14c|ll4cfjl . we obtain the spectral equation. 

We start from equations ()14a|) . Let us introduce new variables A™, i = 1, . . . N, 
associated with the steps of coefficients A and B at singular points, 

A- = {AU - A-)I n {x t ) = (B? - BL x )K n {xi). (15) 

Obviously As and S's may be expressed in terms of A's: 

. A™ 



,)- 



-X ^ny-^j, 



On the other hand this representation reduces the number of unknowns by N and 
guarantees that equations (|14ajl are fulfilled. Note that all A™ and B™ depend only on A" 
with smaller numbers j < i. 

Substituting definition (|T5|) into ()14b|) we obtain the following equation for A, 

A „ = 2K n (x i )I n (x i ) [A?I n ( Xi ) + B?K n ( Xi )] 
l ~ KX[K n ( Xi K( Xi ) - K^x^Uxi)} • l ' 

The expression in square brackets in the denominator is nothing but Wronski determinant 
of the modified Bessel functions I n and K n 

W(x) = K n (x)I' n (x) - K' n (x)I n (x) = -. (18) 

x 



Inserting into (J17J1 formulae (|16j) for A" and Bf we obtain after obvious cancellations: 

t-i 



a n ^X {t\ n \X {j 1 n \X i) 

1 ~~ = ^A 



^M + ^M + EAjgg-lg) 



(19) 



Thus we managed to express every A" in terms of the preceding A™ with j < i. 

As long as we are interested in the spectrum we may simplify calculations by taking 
Aq = 1, Bq = in accordance with (j!4cj) . (The correct value of Aq can be found later 
from the normalization condition (fTH|) .) Given these values we can find A™ and then 
calculate the successive values of A™ A™, £?", i = 1, . . . N. Every new tube brings an 
extra power of (kA) _1 . Thus after i recursions we obtain a polynomial of power i in -r. 

The last thing is to impose the boundary condition (j!4d|) . This determines the spec- 
trum of the system. 

We have already explained that A 7 ^ is an iV-th order polynomial in -A-. However its 
coefficients contain Bessel functions that depend on x\ = kti and one may not be sure 
that the equation A^(k, A) = always has N roots in k,. 1 Still, as we shall see later, the 
roots lie close to each other and variations of Bessel functions are negligible. Hence in a 
system of N coaxial nanotubes each radial energy level is, in general, split into N ones. 

3 The spectral equations 

3.1 The single tube 

Let us analyze the electronic spectrum for a nanotube of radius r t . Taking the wave 
function in form (jHJ) we obtain the equation for radial part: 

2 
x 2 ip'n + xipL — (x 2 + n 2 )ip n = -x 2 6(x — x t ) 4> n , where x t = nr t . (20) 

KA 

The wave function inside and outside the tube looks as follows: 



n{x) = 


= A%I n {x) + B%K n (x) 


for 


X <C Xt] 


(21a) 


n{x) = 


-- A n x I n {x) + B?K n (x) 


for 


X > x t - 


(21b) 



Coefficients A and B must satisfy equations (|T4"j) with N — 1. As follows from the 
continuity condition (j!4a)l we may introduce an auxiliary variable A" such that 

A? = (A n - Al)l n {x t ) = {B? - BX)K n (x t ). (22) 

For the solitary potential well equation (JT§j) is very simple and takes the form 

= 2x t K n (x t )I n (x t ) + B n Kn{xt)] (23) 

KA 

Now we should impose the boundary conditions and demand the wave function to be 
regular at x = and x = oo. This means that 

B% = A\ = 0. (24) 

1 We may suggest, though, that for any k there are N values of Xi such that A^(n, Ai) = 0. 



We have already mentioned that when searching for the spectrum we may assume that 



Aq = 1. Equating expressions (|2*2j) and 



—K n (x t )In{xt) 



for A" we arrive at the spectral equation: 

1 (25a) 



or in terms of the original physical radius: 

2r t 



A 



K n Ur t )I n Ur t ) 



(25b) 



Solution of this equation is greatly simplified by the presence of small parameter 
\/r t <C 1. This justifies the use of asymptotic expressions for the Bessel functions: 



I n (x) 


r~^ 


e x 

\ / 2tcx 


1 


K n (x) 


r^ 


\l^ 


X 



— {n 2 --) + —{n 2 --){n 



2x 



A 



6X^ 



A' 



1 ~)+o(x- 2 ) 



1 + —in 2 --) + -±-(n 2 - -)(n 2 --) + o( X - 2 ) 
2x K A J 8x 2K 4 A A J v ; 



(26a) 
(26b) 



Substituting these into ()25aj) we get with (nr t ) 2 accuracy: 



k\ 



2x 2 



[n 



h 

A' 



(27) 



There must be three roots to this cubic in k equation. Two of them are small (k 2 ~ 
{n 2 — ~)/2r 2 and lie out of the domain where asymptotic formulae (J26j) hold. The third 
one can be found perturbatively. Taking k = k^ + k^ 1 ) with k^ 1 ) <C «(°) = A -1 we get for 
the first approximation: 



k = T 



2r\ 



:n 



1, 

4' 



(28) 



Note that: a) k ~ A -1 is big; b) the correction \{n 2 — \)/2r 2 is small; c) x t ~ r t /X is big 
and we may use the asymptotical form (|26|) of Bessel functions. Thus our approximation 
is consistent. 

According to equation © K determines the electronic energy spectrum: 



E n {p z 



Vz_ 

2m 



H 2 



2mA' 



+ 



h 2 



2mr 2 



n 



mU z 



2h 2 8 



+ 



mr: 



2m 



2^2 



+ 



h n 
2mr 2 



mil 2 

' 2h 2 



(29) 



is the 



The origin of all the four addends is transparent. The negative first term 

binding energy of electron in the potential well. It does not depend on the form of graphite 

layer. The last two terms are the kinetic energy of the electron: the longitudinal part 



r, z , 



2m 



varies continuously whereas the angular component ^^ is quantized. For the tube of 

radius r t = hnm the scale of quantization is -^—>x ~ 1.5 meV . An interesting feature of 

the spectrum is the geometrical term — -^ that depends only on the radius of the tube. 
Because of it electrons are bound to the tube stronger than to the plane. The thinner the 
tube, the stronger the binding. For r t = 5nm it is h 2 /8mr 2 p» 0.38 meV). 



3.2 Two tubes 

A new phenomenon that occurs in multiple nanotubes is the electron tunnelling. In planar 
systems the hopping of electrons between layers leads to splitting of degenerate energy 
levels. Out of the planes the wave functions of electrons fall exponentially and the splitting 
is weak compared to the binding energy. However the latter is large by itself and the effect 
might be noticeable. Now we are going to consider the tunnelling in coaxial nanotubes. 
Of course the presence of curvature terms (oc — ^-j) eliminates the degeneracy from the 
very beginning. Still tunnelling shifts energy levels and enhances their splitting. First we 
shall analyze this effect in the system of two coaxial nanotubes. 

Let the radii of the tubes be r\ and r 2 . The parameter characteristic of tunnelling 
is e~ a = exp (— na) ~ exp (—a/A), where a = r 2 — r% is the interval between the tubes. 
We shall analyze the electronic spectrum of the system in the first approximation with 
respect both to tunnelling and localization of electrons. Thus we will have at our disposal 
two small parameters e~ a and A/ri,2- 

Following the recipe of Section 12. HI we assume that Aq — 1, Bq — and calcu- 
late A^{k : A, n, r 2 ). The boundary conditions require it to be zero. Imposing this 
condition gives the spectral equation we are looking for. At first glance the equation 
A%(k, A, ri, r 2 ) = looks rather clumsy but it simplifies after neglecting terms of the 
order e~ 2a jx\ and higher: 

f 1 - ^ - \ (« 2 - 3 (3 + k) (1 - Xk) + ii (" 2 " j) 2 " '"" = a (30) 

The coefficients of the equation depend on k via x~ 2 and e~ 2a . However this dependence is 
weak. We shall consider them as constants and substitute for k the zeroth approximation 
value ft(°) = A -1 . Then the equation becomes a simple quadratic one. The two solutions 
to it are 




(l-AKK^j^-jj^^jiy-^-.j I-,--,] -He- (31) 
Calculating the energy levels from equation © we get: 



h 2 p 2 



Am 



n 2 




1\ 2 /1 1\ 2 16e- 2 « 



A 1 



.(32) 



Let us discuss the limiting cases. Suppose that the attractive potential be strong 
enough to prevent tunnelling. This corresponds to e~ 2a — > and only r~ 2 terms (i. e. the 
rotational part of kinetic energy) should be taken into account. The two roots of equation 
flUD are 



,-! 1 \ c„„ ia„-2a ^ ( „2 ' 



't-^-TzrA" *) for 16e *K n v U"^' ■ '**" 



That corresponds to the energies 

r, ^ & V \ h 2 (n 2 -\) £ 16e~x /„ l\ 2 /l 1^' 

E n ( Pz ) i, 2 =- - 7TZ7 ^ + i±+ \_ 2 4 ; for -^— «(n 2 - 7 ) (-"- 



2mA 2 2m 2mr 2 2 A 4 V 4 



rf r 2 



(33b) 



Comparing this with (|2*7|) we see that the spectrum consists of two independent series 
coming from the noninteracting tubes of radii r\^- This is exactly what one should expect 
in absence of tunnelling. 

The second case corresponds to large radii and small distance between the tubes. Now 
we neglect x\\ <C e" 2a . Now the solutions become 

(1 - A«)i, 2 = ±exp (-a) for 16e- 2a > (n 2 - i J ( \ - -^ J , (34a) 

and the energies are 

^n(p z )i,2 = -7r-T? lT2exp (-- +^ for -T^>U 



2mA 2 V ^ V \) ) 2m A 4 V 4 



(34b) 



This coincides with the splitting of the ground level that arises in two parallel planes 
due to electron tunnelling between them. Note, however, that because of the factor 
(n 2 — j) dropping down x~ 2 terms is legitimate only for low-lying levels whereas for large 
n the kinetic part of (}3*2"j) always comes into play. When neither of the terms under 
the square root is small the two effects add up. Thus we may conclude that tunnelling 
enhances the splitting of equal n energy levels in coaxial systems. This effect is more 
relevant for small n and becomes less for upper levels. 

Summarizing this Section we notice, that the value of tunnel splitting AE tunne i = 
-^j2&~^ depends on distance between the tubes. It reaches the maximum for a = 2A: 

2fi 2 f) 2 

^KZei = -^- 2 e~ 2 « 0.27—. (35) 

In real systems A is on atomic scales (several angstroms) and a is several nanometers. 
Therefore probably the maximum is inaccessible. In order to increase the role of tunnelling 
one has to make the tubes closer to each other. 

We can also estimate the number n starting from which the effect of tunnelling becomes 
less important and the systems behaves like the set of independent tubes. In order to 
obtain the upper bound let us assume that a = 2A <^ ri j2 . Then inequality ()33b|) turns 
into _ 

O , „ 3 , „ 3 

2 / r\2 /r\2 



">VeW =°' 86 W • (36) 

This is a large value. Certainly in real systems with o>A the tubes become independent 



much earlier. Actually for n this big the kinetic part in (|33bj) exceeds the binding energy. 
Therefore the total energy is positive and the entire model fails. An improvement of this 
estimate requires a better knowledge of A. 

9 



We see that the tunnelling of electrons between coaxial nanotubes enhances the split- 
ting of energy levels with equal n. This effect is stronger expressed for low-lying levels. 
For big n the role of tunnelling falls down and the energy spectrum approaches the super- 
position of those of several noninteracting tubes. The crucial parameter for tunnelling is 
exp (—a/A). In our model A is a function of the energy of electron binding to conducting 
layers. The tunnel interaction enhances as the tubes get closer and the interval between 
them becomes small. An estimate for the tunnelling parameter may be obtained from the 
probability of electron "hopping" between the layers in planar systems. 

4 Physical properties 

4.1 The density of states 

As it is well known the one-particle density of states (DOS) can be expressed in terms of 
the imaginary part of one-electron Green function. In the model of MWCN introduced 
above the quasiparticle spectrum has a form: 

2 

£(p z , n ,k) = p- + AE(n, k) - E F , (37) 

where AE(n, k) denotes the energy related to circumferential modes quantization and 
splitting due to the inter-tube tunnelling. It may be numerically calculated for MWNT 
with any number of conducting cylinders with the help of the formalism developed in Sees. 
2 and 3. Index n enumerates the levels of angular quantization, while index k enumerates 
the energy levels related to electron hopping between adjacent tubes (fine structure of 
energy levels). 

The electron Green function is, therefore, 

G R (p z , n, k; e) = - 2 , (38) 

e-£-AE(n,k)+E F + ±' 

(r is the time of one-electron elastic relaxation). 

The DOS of electron in the multi-wall nanotube is now given by 

i/(e) = --Imf^y - 2 . (39) 

71 J-™ 2vr tte-^-AE(n,k)+E F + £ 

The summation over n should be performed to take into account all the occupied electronic 
states, i.e. from — iV max to iV max , where iV max is determined by the value of electron chem- 
ical potential. A typical value for a realistic nanotube is iV max « 5 [T^. The summation 
over k should involve all energy levels of the fine structure. The numbers of those is equal 
to the number of coaxial tubes in MWNT. After integration over p z and straightforward 
algebra we get: 

= ^ V l (40) 

In order to calculate thermodynamic functions in terms of Green function we follow 
the formalism of Ref. ^3j- Inasmuch as the analytical properties of the Green function 

10 



are the same as in the case of normal 3D metal, we may use the result of Ref. ^] after 
summation over e. In the case of MWNT we have: 



S_ 
V 



2ttT f°° dp 
2ttT 



^Vlm 
9ir Z-w 



n.k 



1 



dG R (p z ,n,k;e) 



-Im 



dp 



a Pz y^ 

9ir 2-^ 



_G R (p z ,n,k;e) 
1 



Oe 



3 "J- 27[ ^i + ^E(n,k)-E F -^ 
and we see that the entropy is proportional to DOS: 



e=0 



(41) 



S 4vr 



Ti/(0). 



(42) 



V 3 

4.2 The tunnel current 

We study now the tunnel current between the nanotube and normal metallic contact. The 

general expression for the tunnel current is 



I(V) = -eIm[K R M]^ =eV , (43) 



where 



/°° dr> f d 3 k 

T^E / 7^3 \T p , k \ 2 G(p z ,n,k;e m + u u )G^(k,e m ). (44) 
:„. -°° n,k J ( 2vr ) 

Here the superscript (0) refers to the normal metal. Integration of the corresponding 
Green function over momenta gives 



d 3 k 

(27rf 



G( >(k,e TO ) = -wri/°>sgn(e ro ) 



(45) 



One can also easily perform the p z integration of the Green function G(p z , n, k; e m + u> v ): 



— L G(p z ,n,k;e m + u} l , 

Z7l 



t7c^msgn(e m + lo u ) 



y/2[E F -AE(n,k)+ie m+l/ } 
The calculation of the sum over e m is straightforward: 



(46) 



Im 



y^ sgn(e m ) sgn(e m + u v ) 
e V- AE ( n > k) + E F + ie m+v 



u:a 



iiOv=eV 



. e + eV e 

tanh — — tanh ■ 



(47) 



ZnJ-ao'" V 2T IT) ^/E F -AE(n,k) + e' 

Finally, for the case of low voltage we can expand in eV and obtain the Ohm's law: 



11 



J(V0 = 2e 2 lV°)(0M0)<|T p , k | 2 ). (48) 

For intermediate voltages the non-linear regime in tunnel conductance may be evaluated 
numerically: 



G tun (V) = dJ ^ = eV ^ <|T p , k | 2 ) ^ f 



dI (V) _ „ 2 „ (0) rzz * nrr i2\V- r _Jl 1 . 

( „ i-oo cosh 2 ^ y/ e + E F - AE(n, k) - eV' 

(49) 
When the temperature is not too low (E F — AE(n max , A; max ) « T « y^z) one can 

substitute the cosh 2 ^ by unity and using the fact that the distance between the levels 
of angular quantization is still much larger than the temperature one gets 

G tun (V) » e 2 N tub u^V^^ <|7p,k| 2 ) V2T - eV. (50) 

In the limit of very low temperatures T <C E F — Ai?(n max , k max ) the cosh" 2 ^ works as 
the delta-function and the tunnelling conductance reproduces the energy dependence of 
density of states: 

(51) 
This result is valid for temperatures less than the fine energy-level splitting T <^ 
^exp(-f)~lX. 

5 Summary 

Let us summarize the results. Using the Kronig-Penney type potential we model the 
MWNT as a cylindrical crystal. Because of the circumferential modes quantization, the 
electron states in the individual tube form a sequence of one-dimensional sub-bands. A 
new phenomenon that occurs in multiple nanotubes is the electron tunnelling between the 
adjacent tubes. The hopping of electrons between the tubes leads to splitting of otherwise 
degenerate energy levels. We have developed the general formalism which allows to eval- 
uate numerically the energy levels of nanotube composed of any number of conducting 
cylinders. We have demonstrated that tunnelling enhances the splitting of energy levels 
with equal n in coaxial systems. This effect is more relevant for small n and becomes 
less for upper levels. The crucial parameter for tunnelling is exp (—a/A) , where a is the 
distance between the tubes and A is the effective width of 5-potential. Naturally, the tun- 
nel interaction enhances as the tubes get closer and the interval between them becomes 
small. 

We applied the Green function formalism and calculated the electron spectrum, the 
one-particle density of states and entropy. Moreover, we demonstrated the possibility of 
direct experimental measurement of DOS structure by means of the tunnelling microscopy. 
The calculated within our approach theoretical dependence of tunnel conductance on 
voltage is presented for realistic set of physical parameters in Fig. 1. In this plot we 
choose the parameters typical for MWNT: U — 1 Ry (pseudopotential of carbon), a = 0.3 
nm, and R — 17 nm. Notice, that in this case the effective width of attractive singular 
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potential A is 0.02 nm confirming our approximation of strong potential U ^$> -^-y. One 
can see the resemblance between the fine structure of peaks of the theoretical density of 
states and experimentally observed peaks of tunnelling conductance |14j . 

An important issue is the possibility of direct experimental studies of DOS by methods 
of tunnelling spectroscopy According to Eq. (|4*9*)l . the tunnelling conductance is the 
convolution of two functions depending on the electron energy This equation results in 
direct proportionality between tunnelling conductance and DOS in the limiting case of 
very low temperature. Namely, the temperature should be much less than the energy scale 
of DOS to be resolved. Therefore the interpretation of results of tunnelling spectroscopy 
needs care. In particular, for resolving the fine structure of DOS due to the intertube 
hopping the temperature should be not higher that few tenth of Kelvin. 

Our theory is based on the analysis of the non-interacting electron system and does not 
explain the appearance of the experimentally observed low-temperature zero bias anomaly 
in MWNT |15| I12j. Nevertheless, it is clear that such dip in the tunnelling conductance 
can be attributed to the effects of electron-electron interaction either in diffusion channel 
(Altshuler-Aronov-Lee effect) [THIUS] or in Cooper channel (superconducting fluctuations) 
[T5] . In the first scenario the width of the dip on the voltage scale turns out to be of the 
order of t^ y {t is elastic scattering time), while in the second one it is of the order of 
7rT. For the case when zero-bias anomaly is large in magnitude, the nonperturbative 
approaches have to be applied [19J. Anyway, experimental data should reflect both free- 
and interacting electron effects in tunnel resistance. This would lead to a complicated 
dependence on voltage as it was really observed in MWNT. 

Another important aspect of direct probing of intrinsic electronic properties of carbon 
nanotubes concerns the experimental problem of making low-resistance contacts to the 
nanotube. Indeed, if the resistance of the contact between a nanotube and metallic 
electrode is much higher than the resistance of the nanotube itself, the low-temperature 
data would refer to the properties of the metal-nanotube contact, rather than to intrinsic 
properties of the nanotube [HJ E] . 

Authors acknowledge financial support from the Program "New Materials" , "Univer- 
sities of Russia" of Russian Ministry of Education and RFBR grant 03-02-16209 , Italian 
projects FIRB "TIN" , COFIN 2004. 
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Figure 1: The theoretical dependence of tunnel conductance on voltage for realistic set 
of physical parameters. 
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